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f^ . Abstract 

O ■ It can be shown that negative energy requirements within the Alcu- 

^J ' bierre spacetime can be greatly reduced when one introduces a lapse func- 

» I t tion into the Einstein tensor. Thereby reducing the negative energy re- 

6J[). quirements of the warp drive spacetime arbitrarily as a function of A(ct, p) . 

With this function new quantum inequality restrictions are investigated 

. rH , in a general form. Finally a pseudo method for controlling a warp bubble 

j^ ■ at a velocity greater than that of light is presented. 

. S : 

1 Introduction 

In recent years the possibility of interstellar travel within a human lifetime has 
become a hurtle for theoretical physics to over come. The discussion involves 
the use of radically transforming the geometry of spacetime to act as a global 
means of providing apparent Faster Than Light (FTL) travel. In 1994 Miguel 
Alcubierre then of Wales University introduced an arbitrary spacetime func- 
tion which provided an apparent means of FTL travel within the frame work of 
General Relativity (GR) Q. The drawbacks of this revolutionary new form of 
propulsion results from the production of causally disconnected spacetime re- 
gions and major violations of the standard energy conditions within GR. Further 
investigations into the warp drive spacetime have suggested that any manipu- 
lation of the spacetime coordinates in regards to FTL travel requires negative 
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energy densities M. Alcubierre's generic model coupled with causally discon- 
nected spacetime regions and large negative energy densities make it hard to 
build a justifiable model of a working warp drive spacetime. However recently it 
has been shown that quantum inequalities can allow for the existence of negative 
energy densities ||, and it has also been found that classical scalar fields can 
generate large fluxes of negative energy in comparison to Quantum Inequality 
(QI) restrictions ||. Coupled with recent astrophysical and cosmological obser- 
vations in support of negative energy in the form of Quintessence |5J we feel that 
there is still new life to be brought into the warp drive spacetime. The once 
distant dream of reaching the stars may become a reality with the aid of the 
warp drive spacetime using latin we can express our motto as Ex Somnium Ad 
Astra (ESAA) which translates; from a dream to the stars[]. The following work 
was made possible as a collaboration of the ESAA motto in order to produce 
a physically justifiable model of the warp drive spacetimeQ The inspiration 
behind this work was to choose a somewhat less arbitrary means of defining the 
Alcubierre Warp Drive spacetime. Specifically we have set out the task of re- 
defining the functions of the warp drive spacetime in order to reduce the overall 
energy requirements as seen in respect to the Einstein tensor. 

1.1 The Alcubierre spacetime 

The Story of the Warp Drive begins with the Alcubierre paper [Q, where we 
have the following metric: 

ds 2 = di 2 - (dx - v s f(r s )dt) 2 - dy 2 - dz 2 . (1) 

We have also modified the metric signature in (Q) to correspond to (1, —1, —1, —1), 
for reasons which will become clear in later sections. The Alcubierre spacetime 
is further defined with the following functions 

v s (t) = ^ (2) 

r s (t) = [{x-x s (t)) 2 +y 2 + z 2 ] 1 l 2 . (3) 

This spacetime model is often referred to a 'top hat' model, by means of a bump 
parameter a, which is seen through 

_ tanh(cr(r s + R)) - tanh(cr(r s - R)) 
mj " 2tanh(a J R) ( ' 

This equation then allows for a "warp bubble" to develop in which a star ship 
may ride. The generalities of the Alcubierre top hat function can be described 
by the given volume expansion 

9 = v s — -—. (5) 
r s dr s 

lr The authors would like to note that the ESAA motto was created by Simon Jenks who 
helped to formalize the early discussions regarding this work. 

2 The collaboration of this work was made possible through the online discussions forum; 
Alcubierre Warp Drive currently available at the fo llowing URL: 
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This is where we receive our first scientific definition of a 'warp drive.' The 
volume elements expand behind the spaceship, while contracting in front of it. 
The draw back of this work occurs with equation 19, where the energy density 
T ab n a nb violates the weak, strong, and dominate energy conditions as seen by 
an Eulerian observer jjj: 

, = r-V„» = ^r» = ^ G » = -^fj(J-) 3 . (6) 
2 The ESAA spacetime 

One can choose to describe the function of the warp drive in relationship to the 
function g(p), where p — r s and g(p) = 1 — f(p). We approach this aspect in 
a manner not to dissimilar to the one chosen by Hiscock S where we impose a 
pseudo metric transformation through the following functions 

.900 = A{pf - [v s (r) g( P )} 2 (7) 

.9oi = .9io = -v„ g{p) (8) 

.9ii = .922 = .933 = -1 (9) 

In order to reduce the energy requirements arbitrarily we consider the following 
lapse function A(ct, p). Where have set A = 1 both at the location of the ship, 
and far from it, but allow it to become large in the warped region following the 
lapse function A(ct,p). Such that the above functions (0-^) have the following 
metric signature in cylindrical coordinates: 

ds 2 = [(A) - v s (r) g(p)) 2 + 2v s g(p)dctdz' - dz' 2 - dr 2 - r 2 d<j) 2 . (10) 

We further define (|lij) with the following coordinate transformation: 

ct 

v s dct (11) 

such that the ship velocity has the following definition v s = 1 = c thereby 
arriving at 

ds 2 = {A(ct, pf - v s {ct) 2 g(p) 2 )dct 2 - 2v s (ct)g{p) cos(6)dctdr 

+ 2v s {ct)g{p) sin(0)p dctdO - dr 2 - p 2 d9 2 - p 2 sm 2 (d)d(f> 2 . (12) 

Thus the exact the solution for the Einstein tensor G 00 for (1T3) is given by 

G ctct = -\v s {ct? (f p 9(pj) sin 2 (0)/(A(ct,p) 2 - v s (ct) 2 g(p) 2 

+ v s {ct) 2 g{p) 2 cos 2 (e)+v s (ct) 2 g(p) 2 sin 2 (0)) 2 (13) 

where we have imposed the trigonometric identity sin 2 (6) + cos (9) = 1. 

We now just briefly discuss how the dramatic energy reductions can occur 
within the warp drive spacetime. First beginning with the Alcubierre metric in 
the ship frame coordinates which is stated so that the function of interest is g(p) 
instead of f(p). Alcubierre originally introduced a gravitational time dilation 
term that doesn't limit ship speed in the goo term of his metric that has been 
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overlooked thus far. Instead of inserting this term starting with the remote 
frame metric as Alcubierre did, we have inserted it into the ship frame metric. 
With these coordinates an exact calculation of the entire stress energy tensor 
can be made fairly easy. The T 00 term is the ship frame energy density term 
debated to be to large. Its solution with the lapse function A(ct 1 p) inserted is 

T oo_ v* c 4 fdg\ 2 {sin(6>)} 2 
4 8ttG \dp) {{A{ct,p)Y 

see section to see how this equation was derived. Letting A become large 
arbitrarily reduces this negative ship frame energy density requirement. Due to 
the non vanishing presence of the other stress-energy tensor terms there will still 
be negative energy in other frames, and thus a Weak Energy Condition (WEC) 
violation, but we really don't know that the "quantum implied" weak energy 
condition is a frame covariant principle even though it has been expressed in 
a frame covariant manner. [T a bU a U b > for all time like U a .] Since the ship 
forms the warp, satisfying that there is very little ship frame negative energy 
may be enough. 

3 Quantum Inequalities 

Pfenning applied |7]] a quantum inequality for a free, massless scalar field to 
the Alcubierre warp even though the Alcubierre spacetime is not the result of a 
free, massless scalar field. Because of this, the results for the restriction on the 
warp shell's thickness are unreliable. It is this quantum inequalities restriction 
on the warp shell's thickness that causes the negative energy's magnitude to be 
so great. Therefore this result is also unreliable. Pfenning also did not include 
the possibility of an A(ct, p) term other than A — 1. Even though the quantum 
inequality may not be reliably applicable to the warp drive spacetime we can go 
back and redo the calculation with a variable A(ct, p) included. The quantum 
inequality is 

to [ + °° (T^U,U V ) h 3 

In order to apply it we must find T^U^U" and To for a Eulerian observer. A 
Eulerian observer is an observer who starts out just inside the warp shell at its 
equator with zero initial velocity with a small sampling time according to the 
ship frame who "samples" the time it takes To for the negative energy region to 
pass him. 

From the above interval we have eq. (lOft 

c 2 = c 2 (,4 2 - v 2 s g 2 )(dt/dT) 2 ~ 2v s gc(dt/dT)U - U 2 (16) 

thus from the definition of the lapse function we have 

f T =A~\ (17) 

Combining these (with a minimum sampling time) results in 

= \u 2 + v s gcA- x U + \v 2 s g 2 A- 2 (18) 
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[U"] = 



U = -v s gcA x 

cA- 1 
-VggcA^ 1 



[U M ] = [cA ]. 
By inserting T^U^U V = A 2 T 00 c 2 into the inequality (|f) we have 

3 



TO 
IT 



+oc A 2 T m h 

-s *dr > 



c 3 327r 2 r 4 



and 



TO 



/? 2 /d ff 



^i 2 p 2 v^p 



where we can make the replacement 



d,T 



T 2 + T,f 



2 i 2 

r + z 



< 



Gh 3 

c 7 r 2 Tg 



(19) 
(20) 
(21) 

(22) 

(23) 

(24) 



with r = const. 

Pfenning then makes an approximation of the geodesic motion of the Eule- 
rian observer. Equation eq. 5.11 (note further references to Pfenning in this 
section refer to In). Since we are using the ship frame with A(ct,p) inserted, 
we must include its effect here as well 



-vg(r)A {t)t 







p 2 = r 2 + v 2 g 2 A- 2 r 2 


such that 








to 


r + °° v 2 Ug\ 2 dr ^ Gh 3 
_ 00 A 2 r 2 + v 2 g 2 r 2 \dp ) r 2 + t§ ~ c 5 t 2 Tq 



We make a definition corresponding to Pfenning's equation 5.15 

?• 



V> 



[vg(r)] 



And according to Pefenning 5.16 the integral should be approximately 

Gh 3A 2 



< 



Aip(r + Aip) " c 5 v 2 tp 2 T$ 



So the inequality becomes 



n Gh A 2 [A(r)} 2 

3 r& 7i2^-4 



rA(r) 



;(r) + 1 



(25) 
(26) 

(27) 
(28) 

(29) 

(30) 



This result is the same as Pfenning eq 5.17 with the exception that A(r) is not 
necessarily 1. Pfenning then asserts that this result must hold for sample times 
that are small compared to the square root inverse of the largest magnitude of 
the Riemann tensor components as calculated for the local observers frame. 
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If for simplicity A is kept approximately constant through the negative en- 
ergy region, then the largest Riemann tensor component for this spacetime is 

\R p 00p \*P 2 (dg/dpf(v/A) 2 (31) 

So in this case, the sampling time must be restricted to 

t = aA/v (32) 

where < a < 1. 

Inserting this into the inequality and looking at the case of large A approx- 
imately constant Aq through the negative energy region leads to 

A < (3/Tt) 1/2 (Gh/2Trc 3 ) 1/2 (v/c)A /a 2 (33) 

Choosing a = 0.10, this approximates to 

A < l0 2 {v/c)l Pl A (34) 

Now we see that letting A become arbitrarily large also arbitrarily thickens the 
minimum warp shell thickness. Therefore the -0.068 solar mass calculation |J7) 
had a reachable shell thickness. All that remains then is to divide the Pfen- 
ning -0.068 solar mass result by an Aq which would allow the thickness chosen, 
and which will lower the energy magnitude even farther by several orders of 
magnitude. 

4 Total Energy Calculations 

The energy requirements for the warp drive spacetime can be calculated from 

E = f T 00 dV (35) 

For T 00 we have 

00 _v^ I dg\ 2 _jf_ 1_ 

32ttG \dp) r 2 + z 2 A i { ' 

Where we assume A = const to be large, one can also see that this is identical 
to eq. (pT[). We now choose 

dV = p 2 sin(6)dpd6d(j) 

so that we can write 



which reduces to 



v s c . I dg\ / 1 \ , 



Eesaa =-j f2G{fp) ur^- (38) 

Thus from section we can show similarly that the energy can be given from: 

°°v 2 c 4 / 1 \ 2 / 1 



E =- L Ma J U )p ^ (39) 
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5 Post relativistic warp drives? 

One of the requirements for a warp drive spacetime which has a velocity greater 
than that of light is the existence of negative energy . However this negative 
energy, exotic matter, quintessence, or however you chose to define it can be 
made to pseudoly appear from negative brane tension || . More specifically for 
our case, an Anti de Sitter (AdS) spacetime would b e a logical choice for such 
an exploration. The reason for this choice is the fact that A acts as a scalar 
field to exploit negative energy densities B . In reference to H we can write an 
arbitrary metric for a single "warped" extra dimension that has 3-D rotational 
invariance in AdS space by 

ds 2 = -a 2 (r, t)dt 2 + b 2 (r, t)dx 2 + c 2 (r, t)dr 2 (40) 

This equation at a first approximation is identical to eq. (|1C|). Now considering 
eq. ( |12|) in reference to eq. (|40|) a five-dimensional AdS warp drive spacetime 
can be given by 

ds 2 = -A bK g(p) 2 dt 2 + A(ct, p)d£ 2 K + (A(ci,p)-v s (ci,p) 2 g{p) 2 ) 

dci 2 + dp 2 + p 2 d9 2 + p sin 2 6d<t> 2 (41) 

where 

< = <id£v)+^ < 42 » 

with k being the curvature scale and L being the length scale of the bulk Ab K . 
Where ct, p = const and cE is a unit metric for a 3 brane. From this the 
Alcubierre warp drive function f(p) expands and contracts within the (3+1) 
slice of the AdS spacetime via a scale factor dependent on time (ci, 0). Since 
the five-dimensional AdS spacetime contains slices of the four-dimensional space 
the time coordinates rescale differently at different points in the extra dimension 
of the de Sitter space. Thus causing the (ct, p) terms to scale differently in the 
5-D space than it does in the 4-D space, therefore a gravitational wave in 4- 
D space can appear to exceed the light because the bulk spacetime acquires 
different velocities for c = 1 when moving through the 'extra' dimension. This 
spacetime would also act to suppress the "need one to make one" paradox of 
pi[ . When the 5-D space is static and there is no gauge field one can have a 
more general solution 

ds 2 = -A(ct, p)dt 2 + L- 2 p 2 dY? K + (A(ct, p) - v s (ct, p) 2 g(p) 2 )dp 2 (43) 

with 

L- 2 = -\nlHk (44) 

which would essentially describe a warp drive in which vq < c. What is inter- 
esting to note about eq. (|dj) is that it can make the 4-D warp drive function 
g(p), appear to be generated from a cosmologic term of order Ag(p) in AdS 
space. Meaning that matter within the warp bubble wouldn't need to become 
tachyonic to solve the horizon [control] problem of the Alcubierre Warp Drive 
as c < wo (i-e. the scalar field acts to create a hyper warp drive without large 
negative energy requirements), it is however noted that inducing the fifth di- 
mension is purely a mathematical trick which may have problems of its own. In 
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the tradition of Alcubierre the warp drive of eq. (^Tj) just begs to be named after 
its familiar counterpart, the "hyperdrive" of science fiction. We also note that 
a similar but more conventional approach was taken by Gonzalez-Diaz pfl|, in 
which a Alcubierre-Hiscock spacetime is embedded in a three-dimensional Mis- 
ner space, which fits into the frame work of classical GR as opposed to this 
discussion. 

5.1 a Broeck extended AdS spacetime? 

Equation (|4l|), shouldn't look to surprising at first hand, let us look at the 
Broeck spacetime |Q : 

ds 2 = -dt 2 + B 2 {p)[dx - v s {t)f{p)dtf + dy 2 + dz 2 \ (45) 

We see that L~ 2 p 2 = B 2 {p), thus M3) can be interpreted as an AdS Broeck 



spacetime, this is because B 2 (p) would naturally arise from primed coordinate 
transformation s |13| , e q.(|ll[ ) is a generic example of this. Thus when comparing 
equation 5 of |l2(] to ([43|), something quite odd occurs to the lapse function 
A{ct, p) when vq > c effectively increasing the apparent amount of negative 
energy within the local brane. Moreover a similar scenario for (J43J) implies a 
coefficient of order A 2 -A^gip) 2 - The resulting total energy of the bulk spacetime 
in comparison to (|38| ) when vo > c is thus 

Notice that the bulk comsomologic term Af,fc acts to flip the negative energy 
our local brane. This can be expected because a brane with negative tension 
can contribute and act positively within a higher dimensional brane B, and 
therefore appear to violate the NEC (Null Energy Condition). This also suggest 
that p = k^ = K+ — K~, i.e. the center of the Broeck bubble contains a 
gravitational kink. 

We now want to discuss a specific boundary condition to illustrate the prop- 
erties of a hyperdrive. To begin this discussion the lapse function will be given 
as 

A(ct, p) = apwL = Jpwl (47) 

thus to arrive at specific boundary condition we replace ([f|) with a function that 
resembles eq. 4 of [pi therefore we have the following piecewise function: 



1 p<R 



A 

2 



fpwL( P )={ (l + [(p-R-^)( P -R+^)]/D) N i?-f< /0 <i?+f 

1 p> R+f 

(48) 
It is also noted that the i?± (A/2) term could be a boundary equal to or greater 
than one corresponding to some value a > 1 in standard units) and that this 
value was chosen arbitrarily, however for this instance we have defined it in 
manner which coincides with Pfenning's analysis throughout this work. N is 
the warp parameter which describes the magnitude of the warp factor: 

N(p) warp = B(p) = l + a (49) 
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Figure 1: The ESAA spacetime represented one dimensionally in Cartesian 
coordinates, with ctpwL — 1 this spacetime appears hyper-dimensional in com- 
parison to the Alcubierre spacetime f(p) = 0. This figure represents the varying 
warp parameter N the minimum value N = B(p) = 1 gives the Alcubierre top 
hat f(p) = 1 (black). The median (gray) and maximum (white) represents an 
energy reduction due to a decreased volume as a factor of the warped region. For 
an arbitrary boundary condition apwh = (1 + [{p — R — ^;){p — R '■+ ^)]/D) N •> 
when R = 100 m the "maximum warp" factor is N — 9.3, where the lapse 
function becomes singular. 



this definition describes how the Alcubierre spacetime differs from our owifl. 

With the boundary condition ( |48| ) the energy requirements for this warp 
drive spacetime can be given from: 

T ap n a n f3 = a 2 PWL T 00 

_(i?+A)2iV G 00 

= -[(R+^ N ][v^/A][ Sl n^}[dg(p)/dp} 2 (50) 



With N = 1 eq. @ becomes 



4 



144G (R + i A) 



2 A 



E bk = ::;„,* : t ^ r / p 4 d P (5i) 



(I 



or 

-^ « 1.7 x 10 38 5 = 1.2 x 10 5 • Msun- (52) 



c* 



The ESAA total energy Eesaa is given with N = 1 since vo = c, this depicts 
an Alcubierre like spacetime from footnote y we see that 



oo 



N = / lim a + 1 = 0+1 = 1. 

However, we now have vq > c, thus a > 1 in standard units, by comparing eq. 
5q) to eq. (plh we see that a = 1, therefore meaning the hyperdrive is defined 



3 Where B(p) is derived from eq. 5 of |12J and the warp parameter N we can also define 
alpha as a = N — 1. 
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through the following parameter 



N = 



1 = 1 + 1 = 2. 



Thus when applying an arbitrary boundary condition ctpwh = (b) N for the 
warped region, we see that N, would reduce the energy requirements for an 
arbitrary total energy calculation such as ( J3q ) . Thus it is seen that under these 
conditions the velocity of the warp bubble is to be implicitly interpreted as 
v = 10 2 c. 

In other words using a Pfenning inspired boundary condition it takes just 
under one earth mass to achieve a velocity 100 times that of light. From the 
ship's frame the apparent velocity is much less from [h3| : 



dx' dx/dt/3 

Utt* ~ 1 - fidx/dt 



(53) 



we see that a fifteen day trip to alpha centauri, would give a local apparent 
velocity of v = 0.01002c, thus compounding the aberration effects of Q. 

5.2 metric patching 

From this it can be seen why the warp drive violates the standard energy con- 
ditions in GR, assuming a negative brane tension one has m: 



T^kpku = -ADg^ aeed k ll k u S n -'(n a ) = A D 



n—p 



E «M 2 



S n - p (n a ) (54) 



from this the NEC becomes T^k^ku < 0. When a killing vector K^ is inserted 
in (pi), for the right hand side we have 



T^kaKK^ = aT^ v = A D 



n—p 



£«) 



AM 2 



'(»°) 



(55) 



which is exactly what we would expect from (Ell) , with the presence of a gauge 
field. In terms of the AdS geometry this is caused by the scale factor (ci,p), 
such that the time lapse function is given by the relation A 2 K\,k = A(ct, p) 2 , i.e. 
converting the brane negative energy into positive energy in the bulk frame, as 
suggested above. 

The total energy in the bulk is positive with r s £ [oo, a(t)] thus with metric 
patching the energy in our three brane can be calculated from: 



QM* = ^v + n nu + Ri w _ R 



-//// 



(56) 



Throughout our analysis of the hyperdrive we have considered a metric patching 
of (3+1) within an embedding in E 4 , while in previous works the warp drive has 
been considered in the form (2+1) embedded in E 3 pQ]. A similar construction 
has also been considered for transversible wormholes M in reference to string 
theory, which conforms to an (n+1) bulk and an ([n-l] + l) brane, again as 
suggest by Everett the warp drive appears as a special case wormhole (with the 
exception that this wormhole appears to be a shortcut through hyperspace). 
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6 Summary 

The goal we set out was to reduce the energy requirements of the warp drive 
spacetime to more physically realistic amounts. Our analysis has shown that 
simply considering a warp drive spacetime with an arbitrary lapse function can 
dramatically lower the energy requirements within the Einstein tensor. We have 
shown that Quantum Inequalities are a poor choice for 'curved spacetimes,' but 
even so minor modifications of QFs can also reduce their allowable energy re- 
strictions. Finally we have derived a mathematical trick which may allow the 
warp drive spacetime to survive within the superluminal range. In closing it 
can be shown that the warp drive spacetime can not be ruled out because of 
'unphysical energy conditions' alone, by simply choosing an arbitrary function 
A(ct, p) the energy reductions become rather dramatic. It is however noted 
that such an arbitrary function is merely a mathematical devise, within a four- 
dimensional spacetime there is no mechanism to set N > 1 as far as we know. 
However, it is somewhat easy to consider N = 2 for a possible "hyperdrive" 
spacetime, with our discussion this would correspond to a five-dimensional AdS 
spacetime (in which a "warp bubble" could be interpreted to move with a ve- 
locity that is 100 times that of light), and without experimental verification of 
the fifth coordinate this may also be considered a clever trick of its own. It is 
therefore noted that the warp drive spacetime appears to have several limita- 
tions within semi-classical GR, e.g. the violation of the Null Energy Condition 
(NEC), therefore it is likely that a more plausible warp drive may reside within 
post relativistic corrections of GR, indeed Ex Somnium Ad Astra. 
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A Selected exact contrivariant solutions of the 
Einstein tensor 



The exact solutions to the ESAA metric ( |12[ ) were produced in Maple VI with 
the aid of GRTensorll and are as follows (and for that reason v s is represented 
by (3 in the following calculations): 



G ctp = \[{AA{ct,p) 2 + 4f3(ct)2g{p) 2 cos 2 (6) - 4(3{ct) 2 g(p) 2 

+ P(ct) 2 g(p)sm 2 (8))(-§- p g(p))p 

+ 4/3(ct) 2 3 (p) sin 2 (9) (f p g{ P )) cos(8)(3(ct)) 

I p (A(ct, p) 2 - /3(ct) 2 g( P ) 2 cos 2 (#) + (3{ct) 2 g(pf sin 2 (0)) 2 (57) 

G» e = \(fiA(ct, pfp (£/3(ct)) + 4A(ct, p) 2 [3(ct) 2 g(p) cos(0) 
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2(3{ct)p (£A(ct, p)) A(ct, p) + 0{ct) A 9{p? 

cos{9) sin 2 {9) [§- p g{p))p 

+ 4/3(ci) 4 5 (p) 3 cos(<9) 3 - i/3(ct) 4 g(pf 

cos(0) + 4/?(ci) 4 ff (p) 3 cos(0) sin 2 (0)) 

(§- p g{p)) sin{9)) {p 2 {A(ct,rho) 2 - f3(ct) 2 g(p) 2 
+ P(ct) 2 g(p) 2 cos 2 )(9)+P(ct) 2 g(p) 2 sin 2 (9)) 2 ) (58) 

B Selected exact covariant solutions of the Ein- 
stein tensor 

Again these are the solutions to the metric (112), with the v s to j3 modifications: 

G,,,, = l(8A(ct,p) (-§- p A(ct,p)) + 4(3(ct) 2 g(p) sin 2 (9) [§- p g{p)f sin 2 (9)p 

8f3(ct) 2 g(p) (-§- p g(p)))p(A(ct,p) 2 - P{ct) 2 g(p) 2 

P(ct) 2 g(p) 2 cos 2 (0) + f3(ct) 2 g(p) 2 sin 2 (9)) (59) 



J pp 



G c 



~\P{ct) 2 {-A P A(ct,pY (£zA(ct lP )) g(p) 2 sin 2 (9) 

A[3{ct) 2 9{P? sin»P (jfrg{pj) A(ct, p) 2 

Ap{ct) 2 g(pf cos 2 (9) sin 2 {9) (fiz 9 {p)) P A(ct,p) 2 

4pA(ct lP ) 2 [3(ct) 2 g(p)3 (fitg{p)) sin 2 (9) 
4p/?(ci) 2 . 9 (p) 4 sin 2 (9)A(ct,p) (j^A{cb,p)) cos 2 (9) 

4g(p) sin 2 (9)p( : § I g{p))A{ct 1 pf 

Af3{ct) 2 g(pf cos 2 {9) sin 2 {9) (fi,g{ P )) pA{ct,p) 2 

ApA{ct,p)P{ct) 2 g{pf (§- p g{p)) (f p A(ct,p)) sin 2 (9) 

3p/3(ci) 4 3 (p) 4 sin 2 (9) (£g(p)) 2 cos 2 (9) 

Af3(ct) 2 g(p) 4 shx 2 (d)A(ct,p)^A(ct,p)) 

4f3(ct) 2 g(p) 3 sin 2 (9) A(ct,p) 2 (l p g(p)) 

4/?(ci) 2 5 (p) 4 sin 4 (9)A(ct, p) (f p A(ct, p)) 

8. 9 (p) 4 cos 2 (9)P(ct) 2 A(ct, p) (-§- p A(ct, p)) 

8. 9 (p) 4 cos 4 (9) fJ(ct) 2 A(ct, p) (-§- p A(ct, p)) 

12 5 (p) 4 cos 2 W/3( C t) 2 A( C t,p) (-§- p A(ct,p)) sin 2 (9) 
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8g(pf cos 2 (9)(3{ctf (-§- p g{p)) A(ct, pf 

3[3(ctfsin\9) (jLg{p)f pA{ct,pfg{pf 

3(3(ctf sin 2 {9) (lg( P )) 2 pg(pf cos 2 (9)A(ct, pf 

+ 4pf3(ct) 2 g(pr (-§- p A(ct,p)y sm 2 (8) 

3p/3(ctfg(pf (f p g(p)) 2 sin 2 (6) 



ApP(ctfg{pf(§- p A(ct,p)) 
+ 3p[3(ctfg(pfsin\9)(jtg(p)) 2 

+ 6pA(ct,pfj3(ctf (-§- p g(p)) 2 g(p) 2 sin 2 (0) 

4pg(p) sin 2 A(ct, pf (f- p A{ct, p)) (j- p g(pj) 
+ 4pf3(ctfg(pf sin 2 {9)A{ct, p) (j- p A{ct, p)) (§- p g(p)) cos 2 (6) 
+ ApP(ctfg{pf sin\9)A{ct 1 p) (g- p A(ct, p)) (-§- p g{p)) 

4pp(ctfg(pf sin(9f (f p A(ct, p)) ' cos 2 (8) 
+ 8f3(ctfg(pf cos\9) (-§- p9 (p)) A(ct, pf 

+ 12f3(ctfg(pf cos 2 (9) (f p g(p)) A(ct, pf sin 2 (9) 

+ Wfg(pf sin» (fab)) A^ pf 

8g(pfcos 2 (9)A(ct 1 pf[§- p A(ct lP )) 
Ag{pf sin 2 A(ct,pf(f p A(ct,pj) 
+ sin 2 (9) (£g(p)) 2 pA(ct,pf 

+ 8g(p) cos 2 (9) (f p g(p))A(ct,pf 

+ 4g(p) sin 2 (9) (-§- p g(p))A(ct,pf) 

I (A(ct, pf - li{cifg{pf + (3{ctfg{pf cos 2 (9) 

+ P(ctfg(pf sin 2 (9)) 2 p (60) 
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